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ABSTRACT 

Recent studies have shown that massive galaxies in the distant universe are surprisingly compact, with typical 
sizes about a factor of three smaller than equally massive galaxies in the nearby universe. It has been suggested 
that these massive galaxies grow into systems resembling nearby galaxies through a series of minor mergers. 
In this model the size growth of galaxies is an inherently stochastic process, and the resulting size-luminosity 
relationship is expected to have considerable environmentally-dependent scatter To test whether minor mergers 
can explain the size growth in massive galaxies, we have closely examined the scatter in the size-luminosity 
relation of nearby elliptical galaxies using a large new database (Nair & Abraham 2010) of accurate visual 
galaxy classifications. We demonstrate that this scatter is much smaller than has been previously assumed, and 
may even be so small as to challenge the plausibility of the merger driven hierarchical models for the formation 
of massive ellipticals. 

Subject headings: galaxies: fundamental parameters — galaxies: photometry — galaxies: structure 



1. INTRODUCTION 

Galaxies exhibit a bewildering variety of shapes and sizes, 
but elliptical galaxies are the simplest of all extragalactic ob- 
jects, and their structural uniformity suggests that they might 
be easier to understand than more complex later type galax- 
ies. Most recent attempts to understand the evolution of ellip- 
tical galaxies have been in the context of hierarchical models 
for structure formation (Toom re|[T977l [White & Rees|[T978t 
Khochfar & Burkert 2 003; Bower et ail.||2006[ |De Lucia et| 



al. 2006 Hopkins et al. 



In these models the for- 



mergers do not form many new stars; instead they simply re- 
arrange the existing stars. It has been argued that such dry 
merging may also be responsible for the observed large size 
growth in massive, compact, elhptical galaxies over the last 10 
Gyrs ( Cimatti et al.'2004'; Daddi et al.'2005'; Longhetti et al.' 
2007; van Dokkum et al. 2008; McGrath et al. 2008; Dam- 
janov et al . 2009, Rettura et al. 2010; Hop kins et al.||2009j 



Naabetal. 2009 . van der Wei et al. 2008 ), though this scenario 
requires considerable fine tu ning to reproduce local scaling 
relations (i Nipoti et al.|2009| . This growth is also expected to 
be dependent on environment, with galaxies in higher density 
environments undergoing more mergers and hence showing a 
larger scatter in their size-luminosity relation. 
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mation of elliptical galaxies is linked to mergers, where the l_ 
galaxy's environment is the central paramet er which deter- 
mines its merger history ( |Mo & WIiite| |1996 ). The predom- 
inance of old stellar populations in nearby massive elliptical 
galaxies is inconsistent with large episodes of recent star for- 
mation, which suggests that the mergers which formed ellipti- 
cal galaxies either occurred very long ago (at redshifts z >2), 
or else that they can best be described as the coales cence of 
pre-existing old stellar populations ('dry mergers ') (Thomas 
eraL][2005|rRenzini|[2006| [Graves et aL ''2009a b). Such dry 



Recent studies ( [Gavazzi et al.|[I996 [van den Bergh[[20(j8 
[Nair et al.|20I0| l have found that the size-luminosity relation 
provides the tightest of all the purely photometric correlations 
used to characterize galaxies. However, the environmental de- 
pendence of this relation is highly debated with some groups 
finding a large scatter, as well as a strong environment de- 
pendent curvature in the size-luminosity relation for elhptical 
galaxies (iGiuric in et al.|I988|[Shen et al.|2003)[Ber nardi et al. 
2007 ; DesrochesTt al.|2007[|von der Linden et al.|20 07 , Hyde 
^Bernardi 2009a'), whil e others do n ot ( |Guo et al.||2009r 
Weinmann et al. 2009; Nair et al.|2010(). The purpose of the 
present paper is to show that the size-luminosity relation of el- 
liptical galaxies is well defined by a fundamental line with no 
environmental dependence. We demonstrate that this scatter 
is much smaller than has been previously assumed, and may 
even be so small as to challenge the plausibility of the merger 
driven hierarchical models for the formation of massive ellip- 
ticals. 

Throughout our paper we assume a flat dark energy- 
dominated cosmology with h = 0.7, — 0.3 and J7a=0.7. 

2. SAMPLE 

2.1. Main Sample 

Our sample of 2861 elliptical galaxies is a subset of the 
14, 034 visually classified bri ght galaxies (model g ' < 16, 
O.Ol < z < 0. 1) presented in |Nair & Abraham l pOTOl ) (NAIO) 
with valid measures of environment (described below) and no 
overlapping companions. The NAIO sample is in turn de- 
rived from the Sloan Digital Sky Survey spectroscopic cat- 
alog (Strauss et al. 2002). Visual classification was carried 
out by one of the authors (RN.) and found to agree with those 
from the RC3 within one Hubble type (on average) for the 
1,793 galaxies in common to both samples. It is important to 
emphasize that, unlike samples derived from automated clas- 
sifications, our sample is expected to contain only minimal 
contamination by SO and Sa galaxies. For this analysis, we 
updated the photometry and spectroscopy to use the SDSS 
Data Release 7 (DR7) derived sizes, luminosities and velocity 
dispersions. 

The sizes of galaxies are parameterized by (a) the radius 
enclosing 90% of the galaxy's light contained within twice 
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the Petxosian radius ( Petrosian|1976 Stoughton et al.|2002| l, 
referred to as Rgo; (b) the corresponding radius containing 
50% of the galaxy's hght, referred to as R^q; and (c) the 
seeing-corrected de Vaucouleurs radius, also known as the 
half light radius, R^ . The de Vaucouleurs radii are corrected 
for the known problem originating in slight errors in SPSS 



sky-subtraction o f large galaxies, as specified by Hyde & 
|Bernardi| ( |2009al l [HE]. 

We explored two estimators of galaxy luminosity. The to- 
tal magnitude of each galaxy was specified by (a) the Pet- 
rosian luminosity and (b) the de Vaucouleurs luminosity, both 
corrected for extinction, k-correction (^Blanton et al. 2005) 



and luminosity evolution ( Hyde & Bemardi„2009a) , and the 
latter also corrected for the known problems of SPSS sky- 
subtraction as prescribed by 'Hyde & Bernardi ( 2009a| l. All 
magnitudes in the present paper are in the AB system. 

It is important to understand the size of the intrinsic scatter 
in our measurements. As will be shown, the tightest relations 
are obtained using Rqq. From simulations, we estimate the 
fractional uncertainties in Rgo to be several times the frac- 
tional uncertainty in the sky level. On SPSS images the sky 
background around galaxies is generally known to about ± 
1%, so the corresponding fractional uncertainties in RgQ es- 
timates are typically 4-5%. Other error terms are added in 
quadrature to this, so that a typical .Rgo measurement carries 
around a 5% uncertainty. The effect o f seeing correcti on on 
i?90 is negligible (see Appendix A in Nair et al.||20lO| l. In 
comparison to the errors on the size measurements, the errors 
on both magnitude estimators are small (< 0.01 magfl. 

Finally, we note that in the plots presented below we iso- 
late subsets of galaxies based on central velocity dispersion. 
These velocity dispersions have been corrected for aperture 
effects as specified by J0rgensen et al. ( 1995). 

2.2. Environment Measures 

To study trends with environment, we use two publicly 
available metrics, an A^th nearest neighbor approach com- 
puted by Baldry et al. (2006), and a group catalog algo- 
rithm by Yang et al. (2007). Baldry et al. (2006) measured 
the environmental density for SPSS galaxies with r' < 18, 
0.01 < z < 0.1 and photometrically selected galaxies with 
surface brightness 18.5</ir,5o<24.0. The density is defined 
as E = N/{TTd%), where djv is the projected comoving dis- 
tance (in Mpc) to the Nth nearest neighbor A best estimate 
density (to account for spectroscopic incompleteness) was ob- 
tained by calculating the average density for A^=4 and A^=5 
with spectroscopically confirmed members only and with the 
entire sample. The mean log S for our sample is -0.32. 

Yang et al. (2007) used an iterative halo-based group finder 
on the SPSS NYU-Value Added Galaxy Catalog (Blanton 
[it al. 2005) for objects with r' < 18, and 0.01<z<0.2 with a' 
redshift confidence Cz>0.7. Tentative group members were 
identified using a modified friends-of-friends algorithm. The 
group members were used to determine the group center, size, 
mass, and velocity dispersion. New group memberships were 
determined iteratively based on the halo properties. The final 
catalog yields additional information identifying the bright- 
est galaxy in the group (BCG), the most massive galaxy in 
the group (both used as proxies for central galaxies), esti- 

^ We ignore the systematic offset in magnitude caused by SDSS under- 
estimating sky for briglit galaxies (r< 14mag). However, if we apply the HB 
correction formula (which is an overestimate) directly to the Petrosian lumi- 
nosity, the trends are similar.) 



mated group mass, group luminosity, and halo mass. We use 
the group occupation number as a proxy for environment. 
We are primarily interested in relative evolution between field 
and cluster galaxies and hence define low-density regions as 
galaxies with N<2 and Baldry log S < -0.32, while high- 
density regions are defined as galaxies in groups with more 
than two members A^>2 and Baldry log E > -0.32 (the mean 
N and log E for our sample). 

3. SIZE-LUMINOSITY RELATION OF ELLIPTICAL GALAXIES 

We determined the best-fitting size-luminosity relations 
for our sample by (a) minimizing the scatter in size (direct 
method) and (b) minimizing the scatter orthogonal to the best 
fit line (orthogonal method). While there is a slight change in 
fit parameters between the two methods, our overall conclu- 
sions are not affected. In all cases we found that Petrosian- 
based quantities yielded considerably tighter relations than 
those obtained with the 'traditional' relations defined using 
de Vaucouluers profile fits. 

Figure 1 (top panel) shows the Petrosian Rqq size- Petrosian 
luminosity relation for the 756 elliptical galaxies that are lo- 
cated in sparse environments (left panel), and the 1315 ellipti- 
cals in dense environments (right panel). The points are keyed 
to central velocity dispersion quartiles where black points are 
the galaxies with the highest velocity dispersion in each panel 
and green points have the lowest velocity dispersion. Orange 
and yellow points are the intermediate quartiles. The black 
line shows the best fit relation for all (^3000) elliptical galax- 
ies. The blue lines/regions and red lines/regions in each panel 
denote the best fit (orthogonal) relation in sparse and dense 
environments respectively, with the shaded region denoting 
the 1(7 uncertainties in the slope (determined by 100 boot- 
straps). The parameters of the best fit orthogonal relation 
are inset in each panel and are summarized in Table 1 . The 
lower panels of Figure 1 show the residuals in size about the 
(direct) size-luminosity relation in small bins of luminosity. 
The parameters of the best fit direct relation are inset in each 
panel. The error bars show the median error in size including 
a 1% error in sky. The dashed contours indicate the 25th to 
75th percentile range while the dotted lines indicate the range 
spanned by the 10th to the 90th percentiles. The figure shows 
the following striking features: (1) The radii and the lumi- 
nosities of elliptical galaxies exhibit a very tight power-law 
relationship (linear in log-log space) over a range of ^ 100 
in luminosity. (2) Within the statistical errors the elliptical 
galaxies in dense and in sparse environments appear to follow 
the same power-law relationship, though the scatter in this re- 
lation is slightly smaller in higher density environments. (3) 
The intrinsic scatter in the relationship is comparable to the 
measurement errors in i?9(£] Thus the size-luminosity dia- 
gram of elliptical galaxies defines a 'fundamental line ' in 
log-log space with negligible intrinsic scatter over two orders 
of magnitude in galaxy luminosity. 

It is interesting to compare our results with those obtained 
by Hyde & Bernardi (2009a) [HB], who found curvature in 
the size-luminosity relationship. The HB sample is a factor 
of twenty larger than ours, but the sample was chosen on the 
basis of automated classifications and is contaminated by SO 
and Sa galaxies. The authors analyzed their sample using con- 
ventional galaxy sizes parameterized by the half-light radius 

^ The measurement errors in Rgo are expected to be larger in higher den- 
sity environments where SDSS suffers from over-estimation in sky back- 
ground levels. This error is not included in the estimated errors. 
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Fig. 1. — Petrosian Size-Luminosity relationsliips for nearby elliptical galaxies. The top row shows the relationships obtained from the 
[N^air & Abraham 7201 O i elliptical galaxy sample where the size of each galaxy is parameterized by the radius enclosing 90% of the galaxy's light 
contained within twice the Petrosian radius, J?go (see text for details). The left-hand panel corresponds to galaxies in the field, while the right hand 
panels correspond to galaxies in dense environments. Symbol colours are keyed to corrected central velocity dispersion of the galaxies, in four 
broad bins. In cyan are error bars for four random points. In each panel the best-fit orthogonal linear model (blue/red lines) is superposed on the 
data, with the parameters for the model inset, and a^y are the scatter in the x and y parameters. The black line shows the best-fit relation 
for all elliptical galaxies. The bottom-row shows the residuals in size from the direct best fit analysis, with the parameters for the model inset. The 
dashed lines show the 25th percentile- 75th percentile range while the dotted lines show the 10th percentile- 90th percentile range. The error bars 
show the median error in size, including 1% error in sky, in small bins of luminosity. 



Re- Figure 2 presents the size-luminosity relationship of our 
sample analyzed using the Petrosian half-light radius (i?5o) 
vs. Petrosian luminosity (top row) and vs. de Vaucouleurs 
luminosity (bottom row) keyed to central velocity dispersion 
quartiles. The color coding is the same as in Figure 1. The 
dashed black line shows the best fit relation from Figure 1. A 
number of features are apparent: (1) The slope, scatter and 
environmental dependence of elliptical galaxies using i?5o is 
nearly the same as Rqq. For elliptical galaxies i?5o is just a 
factor of three smaller than i?go. Thus the ratio Rgo/R^o is 
not sensitive to the internal structure (Sersic index) of ellip- 
tical galaxies. (2) The curvature noted by HB (and clearly 
seen in the lower panels) is not seen when using Rqq or R^n I 
This is because R^, the half light radius, is sensitive to the 
profile shape ( Sersic index) of ga laxies, as has been shown in 
Appendix A of |Nair et al. ( 2010| l. (3) The scatter in the size- 
luminosity relation using Rgo (or R50) is ^ 40% lower than 
that obtained with Re. 

4. DISCUSSION 

The size-luminosity relation of elhptical galaxies is closely 
related to th e Fundamental Plane (Djorgovski & Davis 1987 ; 
Faber et al.||1987; Kormendy & Djorgovski 1989) , which de- 
scribes a rather tight relationship between the size, surface 
brightness and central velocity dispersion of elliptical galax- 

* This is also trae when using the complete HB sample which probes 
higher luminosities. 



ies. The 'tilt' (with respect to canonical relationships pre- 
dicted by the virial theorem) and scatter about the Funda- 
mental Plane are thought to be due to variations in either 
age, metallicity or structural non-homology. It is remarkable 
that the scatter in the Petrosian Rgo size-luminosity relation- 
ship (~'0.090 dex in size for all elliptical galaxies) is tighter 
than the scatter recently reported in the Fu ndamental Plane 
(-0.097 dex in size Re ( |Hyde & Bernardi|[r009b )). How is 
it possible that the size, luminosity and velocity dispersion 
information embedded in the fundamental plane formulation 
does not provide a more accurate fit to the observations than 
does the fundamental line defined by the size and luminos- 
ity data alone? Recent work on nearby Virgo cluster galaxies 
( Kormendy et al. 2009| l makes a strong case for the conclusion 
that the internal structure of early-type galaxies depends on 
their minor merger history. If true, one might perhaps expect 
the dispersion about the fundamental plane to depend on the 
galaxy radius that is chosen to define the plane. Figure 2 sug- 
gests that at least some of the tilt, scatter and curvature in the 
conventional fundamental plane may be due to structural (or 
kinematic) non-homology in the elliptical galaxy population 
with Re, and that this may disappear if the fundamental plane 
is defined using a metric of size which is not as sensitive to the 
galaxy profile shape. Note that this suggests that the record of 
the galaxy's merger history may be most strongly imprinted 
in a galaxy's profile, rather than in its overall size, which is 
what is being probed by our investigation. 
It is interesting to consider if the small scatter in the size- 
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Fig. 2. — Decrease in scatter and curvature of the Size-Luminosity Relation obtained by using Petrosian R^o {B.50) sizes instead 
of conventional sizes. The top row shows the relationships obtained using the radii enclosing 50% of the Petrosian flux, Rao vs. Petrosian 
luminosity while the bottom-row shows the relationships obtained using de Vaucouleurs half-light radii, Re, vs. de Vaucouleurs luminosity. The 
left-hand panels correspond to galaxies in the field, while the right hand panels correspond to galaxies in dense environments, such as groups. 
Symbol colors are keyed to corrected central velocity dispersions of the galaxies, in four broad bins. In each panel the best-fit orthogonal linear 
model is superposed on the data, with the parameters for the model inset, o-^^ and a^y are the direct scatter in the x and y parameters. The 
solid black line shows the best-fit relation for all elliptical galaxies. The dashed black line shows the best-fit relation from Figure 1 . The green lines 
indicate the predicted slopes from dry major merging assuming various orbital configurations from Boylan-Kolchin et al. i2006l The steeper line is 
for bound radial orbits, while the shallower line is for the bound orbit with the largest pericentric distance iboyian-Koicnm et ai. 2006 1. 



luminosity relationship for elliptical galaxies poses a chal- 
lenge for theory. The merger theory of elliptical galaxies 
needs to account for both the tightness of the fundamental 
plane (and the tighter size-luminosity relation) and the growth 
of compact elliptical galaxies which account fo r up to ~50% 
of the galax y population at redshifts of 2-3 ( |van Dokkum| 
et al. |2008[ ). We compare our result s to predictions from 
published simulations of dry mergers (Boylan-Kolchin et al. 
[2006 ). Figure 2 shows the predicted size-luminosity relation- 
ships in green for 1:1 mass dry mergers with various orbits 
which preserve a fundamental plane (although not necessar- 
ily with the observed slope). The size calculated by |Boylai>| 
Kolchin et al. ( 2006| ) is not a profile radius but the circular 



aperture enclosing half of the projected stellar mass, i.e. it is 
more comparable to R^q. The steepest relation (slope = 1.2) 
is for the most radial orbit while the shallower relation (slope 
= 0.6) has the largest pericentric separation between the merg- 
ing pair. While the curvature in i?e can be explained by radial 
dry merging, the lack of curvature in Rqq and i?5o further sug- 
gests that Re is more influenced by profile changes induced 
by mergers. The observed scatter (0.09 dex, see Table 1) in 
the Rgo (i?5o)-luminosity relationship seems much too small 
to be consistent with the predictions of dry merger simula- 
tions. In fact, the observed scatter is consistent with the scat- 
ter due solely to orientation effects predicte d by simulations 



suming the galaxies are optically thin). More worryingly, dry 
minor mergers are expected to introduce a similar and possi- 
bly larger scatter in size as do dry major merger s for a similar 
increase in mass (or luminosity) (Nipoti et al. |2009 i. Vari- 
able gas fracti ons of the progenitors c an introduce a further 
scatter in size ('Robertson et al.||20061l. It thus seems highly 



(0.1 d ex with 10"' viewing angles, see Boylan-Kolchin et al. 
2006| l. (Note that the scatter in luminosity is negligible as- 



unlikely that major or minor mergers (either dry or wet) could 
be growing elliptical galaxies while preserving the slope and 
small scatter observed in the size-luminosity relation in both 
low and high density environments. 

The current prevailing theory of size growth of distant 
galaxies can be described as the development of an out er en- 
velope whic h grows about a central dense 'red nugget' (Dam-j 
janov et al.|2O09 ), whose stellar density remains unchanged as 
the galaxy inflates (Hopkins et al. 2009; Bezanson et al. 2009[ 
|Damjanov et al. 2009). The most straightforward expectation 
( [Bezanson et al ||2009| l based on the virial theorem suggests 
that a minor merger will grow the size of a galaxy in direct 
proportion to the additional mass added to the total system by 
the galaxy being absorbed. There will be considerable scatter 
in the growth depending on the relative orbital configurations 
of the galaxies (Boylan-Kolchin et al. 2006). If minor mergers 
are driving this process, not only must the relative masses and 
orbital configurations of the merging systems be fine tuned in 
order to grow galaxies while maintaini ng the slope of the size- 
luminosity relationship at all redshifts ( Nipoti et al. 2009| l, but 
the outcome must also preserve a negligible scatter about the 
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fundamental line over two orders of magnitude in galaxy lu- 
minosi ty. Assuming a sca tter of 0.07 dex in the size-mass re- 
lation, Nipoti et al. ( 2009) find their simulations rule out size 
growth larger than a factor of 1.9 by dry mergers. 

In summary, it is concluded that the elliptical galaxy size- 
luminosity scaling relation determined using Petrosian i?9o(or 
i?5o or Rp) has a much smaller scatter than the same relation 
determined using the half-light radius R^. The relation be- 
tween the Petrosian luminosities and Petrosian-based radii of 
elliptical galaxies is a simple power law: Rqo cx L^-^. The 
dispersion about this 'fundamental line' in log-log space is 
found to be only 0.09 dex in size (0.36 mag in luminosity) 
in the local universe, smaller than that of the fundamental 
plane defined using the effective radius R^. The 'fundamen- 
tal line' appears to be driven mainly by luminosity (or mass) 



and seems to be independent of environment with no curva- 
ture at higher luminosities. The observation that the structural 
properties of elliptical galaxies is both simple and indepen- 
dent of environment suggests that the theory of hierarchical 
growth of elliptical galaxies via mergers is not understood. 
Is our paradigm for galaxy formation merely cracked or is it 
broken? 
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TABLE 1 
Size-Luminosity- Relationships 



Class 


N 


Slope 


Intercept 


Dispersion 


- Orthogonal Fit - 
Dispersion R Dispersion L 






[log(L0)] 


riog(i0)i 

Llog(kpc) J 


[log(i0)] 


[log(kpc)] 


[log(kpc)] 












■'-'petrostan 






All 
Sparse 
Rich 


2861 
756 
1313 


0.596 ± 0.006 
0.583 ±0.011 
0.605 ± 0.013 


-5.256 ±0.062 
-5.123 ±0.117 
-5.349 ±0.138 


0.077 ± 0.001 
0.085 ± 0.003 
0.074 ± 0.003 


0.090 ± 0.002 
0.099 ± 0.003 
0.086 ± 0.003 


0.151 ± 0.003 
0.169 ±0.006 
0.143 ±0.006 








Elliptical galaxies R50 vs. 


^petrosian 






All 
Sparse 
Rich 


2861 
756 
1313 


0.596 ± 0.006 
0.580 ±0.013 
0.616 ±0.013 


-5.746 ± 0.062 
-5.576 ±0.132 
-5.950 ±0.136 


0.080 ± 0.001 
0.088 ± 0.003 
0.075 ± 0.003 


0.093 ± 0.002 
0.102 ±0.003 
0.088 ± 0.004 


0.156 ± 0.003 
0.176 ±0.006 
0.144 ±0.006 








Elliptical galaxies Rp vs. 


-^petrosian 






All 
Sparse 
Rich 


2861 
756 
1313 


0.6114 ±0.007 
0.598 ±0.014 
0.630 ± 0.014 


-5.515 ±0.070 
-5.376 ±0.146 
-5.711 ±0.149 


0.088 ± 0.001 
0.097 ± 0.003 
0.083 ± 0.003 


0.103 ±0.002 
0.113 ±0.004 
0.098 ± 0.004 


0.169 ±0.004 
0.189 ±0.007 
0.156 ±0.008 








Elliptical galaxies Re 


vs. L^eV 






All 
Sparse 
Rich 


2861 
756 
1313 


0.739 ±0.013 
0.695 ± 0.026 
0.815 ±0.023 


-6.992 ±0.140 
-6.538 ±0.265 
-7.792 ±0.239 


0.135 ±0.003 
0.148 ±0.006 
0.123 ±0.005 


0.168 ±0.004 
0.180 ±0.008 
0.158 ± 0.006 


0.227 ± 0.006 
0.259 ±0.011 
0.194 ±0.011 



